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We examine the speed and fidelity of several protocols for state or single excitation transfer in 
finite spin chains subject to diagonal and off-diagonal disorder. We find that, for a given chain 
length and maximal achievable inter-spin exchange {XY) coupling strength, the optimal static spin- 
coupling protocol, implementing the fastest state transfer between the two ends of the chain, is 
more susceptible to off-diagonal [XY coupling) disorder, as compared to a much slower but robust 
adiabatic transfer protocol with time-dependent coupling strengths. 
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I. INTRODUCTION 

Faithful transfer of quantum states between physical 
qubits of an integrated quantum register is one of the 
important prerequisites for scalable quantum computa- 
tion. Typically, qubit-qubit interactions are short range 
and implementing quantum logic gates between qubits 
located at distant sub-registers would involve intercon- 
necting them via quantum channels, or wires which 
may consist of arrays of coupled quantum dots [3, [1] or 
superconducting qubits 0-16 1, atoms in optical lattices 
0,01 J or other realizations of spin chains. 

Quantum channels of permanently coupled spins would 
require no dynamical manipulations during the state 
transfer, but might be susceptible to noise and imper- 
fections. Conversely, dynamically manipulated networks 
can be more robust with respect to certain kinds of dis- 
order, but are more involved requiring time-dependent 
external control. Here we re-consider critically several 
protocols for achieving efficient and dependable — ideally 
perfect — state transfer in disordered spin chains subject 
to physically constrained maximal achievable inter-spin 
coupling rate. The present work is an extension of our 
earlier studies 0, Q on perfect states transfer to more 
realistic scenarios with the aim of quantifying and neu- 
tralizing the influence of static (or slowly changing) noise 
inevitably present in any imperfect physical realization of 
the spin chain resulting in diagonal and off-diagonal dis- 
order. 

After outlining the model, we examine the speed and 
reliability of several state transfer protocols first for ideal 
and then for noisy spin chains, followed by conclusions. 



II. THE MODEL 

The Hamiltonian for a spin chain of length N has a 
general form ^ 
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where ctJ'^'^ are the Pauli spin operators at position j, 
hj determines the energy separation between the spin-up 
and spin-down states playing the role of the local "mag- 
netic field" , and Jj is the nearest-neighbor spin-spin in- 
teraction which can be static or time-dependent. From 
now on we set the anisotropy parameter A = 0; Eq. ([T]) 
reduces then to the Hamiltonian of the XX model, which 
is isomorphic to the Hubbard Hamiltonian for spinless 
fermions or hard-core bosons @, 
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where (cij) is the particle creation (annihilation) op- 
erator at site j with energy hj and Jj now plays the role 
of tunnel couping between adjacent sites j and j + 1. 

Our objective here is to transfer an arbitrary single 
qubit state = a |0) -I- /? |1) between the two ends of 
the spin chain. To that end, we assume that all the spins 
can be prepared in the "ground" state = \0)j and at 
a certain initial time iin = the first site of the chain is 
initialized to \ip)i. Ideal transfer would imply that at a 
well-defined final time tout the last site of the chain is in 
state IV') AT J up to a certain relative phase factor between 
the amplitudes of states |0)Ar and |l)Ar (see below). 

Since the Hamiltonian ([!]) [or ([2])] preserves the number 
of spin- [or particle-] excitations, we need to consider only 
the zero |0) = O^Li |0)j and single excitation \j) = 
|0) [a] |0)] subspace of the total Hilbert space. Then 
the system initially in state IvPin) = Oi\0) + evolves 
in time as |^'(t)) = (7(0 |*in) = a |0)+/3Ef=i W U), 
where U{t) — Texp [-^ j^H{t')dt'^^ is the (time-ordered, 
T) evolution operator. Apparently, only the states in the 
single excitation sub-space { evolve in time with the 
corresponding amplitudes Aj{t) = {j \ U{t) while the 
vacuum (or ground) state |0) remains unchanged. Thus 
perfect state transfer would be achieved for the amplitude 
|^jv(^out)| = 1, provided its phase = arg(Ajv) is fixed 
and known, (j) = (j)o, and therefore can be amended. 

We may quantify the performance of the scheme by 
the transfer fidelity = {ipl pM IV") where pN = 
Tr^( ) = {l-\!i?\AN?) |0)(0| +mMN? |1)(1| + 
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af3*A% |0)(1| + a*l3AN |1)(0| is the reduced density op- 
erator for the A^th site of the chain ^lOj. We then have 
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while the mean transfer fidelity F, obtained by averaging 
over all possible \^p) and after compensating for 0o, 
is given by 
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Thus, for the amplitude I^Arj = 1 but completely random 
phase (j), the fidelity is equal to the classical value of = 
2/3, while for \A]y\ = we have F = 1/2 corresponding 
to a random guess of the qubit state |0) or |1). 



III. STATE TRANSFER PROTOCOLS 

The state or excitation transfer in a spin chain de- 
scribed by Hamiltonian ([1]) [or ([2])] is mediated by the 
nearest- neighbor couplings Jj. Clearly, in any practical 
realization of the spin chain there will be some upper 
limit for achievable coupling strength, Jmax = max{Jj}, 
determined by physical of technological constraints. On 
a fundamental level, this follows from the fact that the 
energy of the system is bounded, which, in turn, limits 
the speed of the state transfer, tout ^ A^/>/max [111- 

A. Noiseless spin chains 

Let us first recall the key facts pertaining to an ide- 
alized spin chain with no disorder. We assume uniform 
on-site energies hj := V j G [Ij Af]j while the individ- 
ual couplings Jj can be freely controlled, subject to the 
constraint Jj < Jmax- 

(a) Perhaps conceptually the most straightforward 
approach to the state transfer between the two ends of 
the chain is to apply a sequence of swap operations im- 
plemented by TT-pulses between the pairs of neighboring 
sites. To that end, with all the couplings Jj set initially 
to zero, we switch on Ji for time ti = 7r/(2Ji), then J2 
for time ^2 = t^/{2J2), etc. till reaching the A^th site. At 
the end of each step, the corresponding state amplitude 
is Aj{tj^i) = —ism{Jj-itj-i) Aj-i(tj-2) = {—iy~^ for 
j ^ 2, . . . ,N. If all the couplings' strengths can be pulsed 
to the maximal possible Jmax, and there are TV — 1 steps, 
the total transfer time is tout = {N — l)7r/(2Jmax) — 
(7r/2)(A''/Jmax) (TV > 1) with the final state ampli- 
tude AN{tout) = i.e., |Ajv(Wt)| = 1 and 
<j)o = (-7r/2)(TV - 1) (mod 27r) . 

(b) We next consider a spin chain with static cou- 
plings Jj arranged in an appropriate way facilitating 
the perfect state (or excitation) transfer. By "static" 
we mean that during the transfer the coupling strengths 
are fixed, but to initiate (at time tin) and to terminate 
(at time tout) the transfer process at least Ji and Jn-i 
should be quickly switched on and off, respectively. (Al- 
ternatively the state initialization of the first site at tin 



and state retrieval from the last site at tout should be ac- 
complished very fast, on a time-scale short compared to 
Ji~A f-i -) Among the many [l[ — in fact, infinitely many 
— possible static protocols for perfect state transfer, 
we focus here on the one proposed in 0, [l3| , and much 
earlier [l3| and in a different context (that of popula- 
tion transfer in laser-driven multilevel atomic or molec- 
ular systems [isjl). which was shown to be the optimal 
one [il| in terms of the transfer time. In this so-called 
spin-coupling protocol, the coupling constants are ar- 
ranged according to Jj = Jq^J{N — which makes 
the s yste m formally analogous to a spin-j7 in a magnetic 
field [lal- This leads to the equidistant energy spectrum 
Afe = 2 Jofe - Jo(A^ + 1) with fc = 1, 2, . . . , AT, and conse- 
quently perfectly periodic oscillations of the single exci- 
tation between the two ends of the chain, according to 
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Thus, at time tout = t^/{2Jq) the amplitude of the final 
state is AN{t) = [-i sin(Jotout)]^"^ = (-«)^"^- Note 
that the strongest coupling is in the center of the chain: 

at j = A^/2 for A^ even, Jn/2 = 5 Jo-^ = Jmax; or at j = 

{N ± l)/2 for N odd, J(W±l)/2 = i JoVTV^ ~ Jmax 

(A^ 3> 1)- Hence, the transfer time expressed through 
is given by tout = (7''/4)(A^/ Jmax), which is twice 
shorter than that for the sequential swap protocol. 

(c) The last protocol that wc consider here is the 
adiabatic state or excitation transfer between the two 
ends of the spin chain using slowly varying couplings Jj 
[3, [l3| - This is analogous to the stimulated Raman adia- 
iDatic passage (STIRAP) techniques extended to mul- 
tilevel atomic or molecular systems [1^. Assume that A^ 
is odd and the individual couplings Jj can be selectively 
and independently manipulated. In the single excitation 
subspace, the Hamiltonian ([T|) [or ^] has an eigenstate 

|*(°)) =-i=[J2 J4 - - - JjV-1 |1) + (-1) Jl J4 - - - Jn-1 |3)+ 
vA/o 



+ (-l)-^JlJ3--- JAr-2|Ar)], 
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with eigenvalue A*^''^ = 0, which is conventionally called 
coherent population trapping (or dark) state [isl . [l9| . 
Thus the amplitude of initial state Ai is proportional 
to the product of all the even-numbered couplings, while 
the amplitude of final state is given by the prod- 
uct of all the odd-numbered couplings, divided by the 
normalization parameter Afo = (J2 J4 - - - Jn-iY + - . . + 
( Jl J3 . . . JAr_2)^. Therefore, if all the even-numbered 
couplings are switched on first, the zero-energy state 
(m would coincide with the initial state |1). This is 
then followed by adiabatically switching-on of all the 
odd-numbered couplings, while the even-numbered cou- 
plings are switched-off, which will result in state ^ 
to rotate towards the final state \N). Assuming that 
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FIG. 1: Dynamics of single excitation transfer in a spin chain of length A'^ = 25 for (a) sequential SWAP, (b) spin-coupling, 
and (c) adiabatic protocols. Top panel corresponds to noiseless spin chains, (Jh ~ o-j = 0, yielding complete excitation transfer 
|Ajv(tout)P = 1 for all cases (a), (b) and (c). The graph above (c) shows the time dependence of couplings Jovcn and Jodd 
normalized to Jmax [cf. Eq. ([S}]. Bottom panel illustrates the results for noisy chains with oh = = 0.15Jmax averaged over 
1000 independent realizations, leading to (|j4jv(tout) P) — 0.2, 0.42 and 0.96 for (a), (b) and (c), respectively. Time is measured 
in units of J^ix a^nd the evolution terminates at the corresponding iout. 



these two families of couplings are described by com- 
mon shape functions, J2, Ja, ■ ■ . , Jn-i — Jcvcn{i) and 
Ji, J3, . . . , Jn-2 — Jodd{t), the amplitudes of the initial 
and the final states are given by 



[^cvcn (^) 



J 



J 



with J\fo{t) = Ef=o[^odd(i)]'"[^ovc„(t)]'(-^-"). Thus, 
complete state or excitation transfer between the two 
ends of the chain can be achieved by applying first the 
^cvcn couplings and then the Jodd couplings, the two sets 
of couplings partially overlapping in time. At time tout, 
when Jodd (tout) > Jcvcn(tout) — 0, the amplitude of the 
final state is AAr(iout) = (-1)'^, i-e., \AN{tout)\ = 1 
and 00 = {—^){^ — l)/2 (mod 2tt). Of course the adi- 
abatic following of the zero-energy eigenstate Q holds 
true if, during the transfer process, the non-adiabatic 
transitions out of j^^"^) are negligible, which requires 
that the rate of change of the coupling strengths be 
small compared to energy separation between j^*^*^)) and 
all the other eigenstates. We can estimate the energy 
separation between the eigenstates in the vicinity of 
maximal overlap between the even and odd couplings, 
Jcvon — Jodd — J- The energy spectrum of the chain 
with homogeneous coupling, J,^^ = J W j G [l,iV], is 
Afe = -2Jcos[kn/{N + 1)] [1 [|. The eigenstate with 
zero energy A^^-' is the one with k = {N + 1)/2 = fco, and 
the nearest eigenstates with indices k = ko ± 1 have en- 



ergies Xko±i = ±2Jsin[7r/(A^ + 1)] ~ ±2Jtt/N {N > 1). 
With J < Jmax, the excitation transfer time tout, being 
roughly equal to the couplings' switching time, should 
then satisfy the condition tout >• N/{2nJ^ax)- 

To summarize the results for noiseless spin chains, the 
transfer time tout for all three protocols scales with the 
number of sites N and the maximal inter-site coupling 
Jmax as (iV/Jmax)- The fastcst is the spin-coupling pro- 
tocol with tout = (7r/4)(iV/Jmax)- It is followed by the se- 
quential SWAP protocol, for which tout — (7r/2)(7V/Jmax)- 
Finally the slowest is the adiabatic protocol tout — 
C(A^/Jmax), with C ~ 8 being a safe estimate [l^l for 
smooth coupling functions that we use: 
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1 



l±erf 



t — 2^out i 2at 



V2, 



(5) 



with at — I tout- Note that the phase of the final 
state amplitude for all three protocols is given by (po — 
(— 7r/2)(7V— 1) (mod 27r), which should be compensated 
for after the transfer. The single excitation transfer for 
all three protocols is illustrated in the top panel of Fig. [1] 



B. Disordered chains 

Employing numerical simulations, we now examine ro- 
bustness of the above described state transfer protocols 
in spin chains with varying degree of disorder. We note 
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FIG. 2: Averaged (over 1000 realizations) fidelity (F) in noisy spin chains for (a) sequential SWAP, (b) spin-coupling, and (c) 
adiabatic protocols. Top panel shows the dependence of {F) on the diagonal disorder ah with aj = Q (upper plots), and on 
the off-diagonal disorder crj with Oh = Q (lower plots), for the chains of lengths A'' = 15, 25, 51. Bottom panel shows {F) versus 
both ah with a,j in spin chains with N = 25. 



that for the spin- coupling scheme, related analysis has 



been performed in [21 



The physical origin of disorder may be two-fold: (i) 
fabrication imperfections of the particular system realiz- 
ing the spin chain, and (ii) noise of the external controls, 
which is assumed to vary slowly enough on the time-scale 
of state transfer tout, as is typically the case in most 
experimental situations pertainin g to coupled quantum 
dots [l^l , superconducting qubits (23l | or atoms [23| . We 
shall distinguish diagonal and off-diagonal disorder. The 
diagonal disorder corresponds to random on-site energies, 
or equivalently the local magnetic fields hj , normally dis- 
tributed around {hj) = with variance a\ (without loss 
of generality, we assume that the energies of the first 
and the last sites of the chain are exempt from disor- 
der, /ii = /ijv = 0, since otherwise the state would 
dephase even before and after the transfer). The off- 
diagonal disorder introduces randomness in the inter-site 
coupling strengths Jj — Jj{l + 5Jj) where 5Jj are nor- 
mally distributed around {SJj) = with variance ctj. 
Consistently with the above description, we will treat 
the disorder as static during each realization of the nu- 
merical experiment for the particular protocol, but com- 
pletely uncorrelated between different realizations. The 
results presented below are obtained by averaging over 
many (typically 1000) independent realizations. 

Figure [T] compares the single excitation transfer for the 
three protocols (a), (b) and (c) in ideal and disordered 



spin chains of length N = 25. In the noiseless chain we 
have perfect transfer |^Ar(tout)P = 1, while in the pres- 
ence of diagonal and off-diagonal disorder characterized 
by standard deviations Ch = crj = 0.15Jniax, the aver- 
aged transfer probabilities are reduced to (|AAr(iout)P) — 
0.2, 0.42 and 0.96 for the cases of (a), (b) and (c), re- 
spectively. Thus, among the three transfer protocols, the 
sequential swap scheme is the most susceptible to noise, 
especially to the off-diagonal disorder which leads to de- 
viations of the subsequent pulse areas from the required 
value of tt; in this particular example, the off-diagonal 
disorder alone is responsible for at least 70% reduction 
of the transfer probability. The spin-coupling scheme is 
somewhat more robust with respect to noise, with both 
diagonal and off-diagonal disorder comparably contribut- 
ing to the reduction of the transfer probability by about 
20% and 40%, respectively. Finally, the adiabatic trans- 
fer scheme is very tolerant to noise, as far as the transfer 
probability is concerned, but is quite slow; in fact it can 
tolerate even more disorder at the expense of slowing it 
further down (equivalent to increasing C). 

The probability of excitation transfer alone is not 
enough to fully characterize the state transfer, since, e.g., 
|^Ar(iout)P = 1 but completely random phase (j) amounts 
to classical information transfer only, and the resulting 
fidelity for quantum state transfer is merely F — 0.66. 
We will therefore quantify the performance of the sys- 
tem subject to varying level of noise using fidelity (F) 
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of Eq. ([3]) averaged over many independent realizations 
of protocols (a), (b) and (c). Figure [5] summarizes the 
results of our numerical simulations for the chains of 
lengths N — 15, 25 and 51. Unsurprisingly, the longer 
the chain the lower the fidelity of the state transfer. We 
find that, for the same values of diagonal an and/or off- 
diagonal a J disorder, the sequential swap scheme yields 
lower fidelity than the spin-coupling scheme. Moreover, 
both schemes are somewhat more susceptible to the off- 
diagonal disorder. The behavior of the fidelity for the adi- 
abatic transfer scheme is, however, profoundly different: 
it is very robust with respect to the off-diagonal disorder 
a J , but much more sensitive to the diagonal disorder ah ■ 
already for an f^ 0.1 the fidelity (F) ~ 0.66 (but then 
decreases slowly with increasing ah)- This is despite the 
fact that the transfer probability (|^Ar(tout)P) remains 
above 0.9 up to ah,j ^ 0.28, i.e., the excitation transfer 
is very efficient up to large values of both diagonal and 
off-diagonal disorder. The adiabatic transfer protocol is 
so sensitive to diagonal disorder because it is slow: during 
the long transfer time iout even little noise in the on-site 
energies ah accumulates to large random phase spread 
over ant out- 

IV. CONCLUSIONS 

We have critically examined the state and excitation 
transfer in disordered spin chains using the sequential 
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